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Abstract 

The  aim  of  this  work  is  to  investigate  the  limit  performance  of  molten  carbonate  fuel  cells  due  to  gas  diffusion  phenomena  in  the  porous 
electrodes  when  high  reactant  utilisation  factors  are  used.  Modelling  and  experimental  activities  are  presented.  An  electrode  kinetic 
micromodel  has  been  performed  accounting  for  the  effects  of  the  diffusion  limits  and  has  been  integrated  into  a  macromodel  in  order  to 
optimise  cell  performance  and  operation.  Specific  tests  have  been  carried  out  to  identify  the  dependence  of  anodic  as  well  as  cathodic 
limiting  current  density  on  operating  parameters  using  planar  square  cells  produced  by  Ansaldo  Ricerche  (ARI).  In  turn,  limiting  current 
measurements  have  been  used  in  order  to  identify  overall  reactant  mass  transfer  coefficients  in  the  complex  structure  of  the  electrodes.  A 
good  agreement  between  simulated  and  experimental  results  is  presented  and  discussed.  ©  2001  Elsevier  Science  B.V.  All  rights  reserved. 
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1.  Introduction 

Today  one  of  the  most  significant  sources  of  pollution  is 
coming  from  energy/power  systems.  In  this  scenario  molten 
carbonate  fuel  cells  (MCFCs)  are  a  promising  solution  for 
the  production  of  clean  energy  with  high  efficiency  and 
environmental  respect,  especially  for  remote  users.  In  par¬ 
ticular,  MCFC  technology  is  now  at  the  stage  of  scaling-up 
to  commercialisation  and  many  developers  world-wide  have 
shown  significant  progress. 

In  order  to  guide  the  development  process  in  terms  of 
optimisation  of  operating  conditions  and  technological  solu¬ 
tions,  a  simulation  activity  has  been  carried  out  and  pre¬ 
sented  in  a  previous  paper  [1];  an  MCFC  macroscopic  model 
has  been  set-up  where  a  micromodel  of  the  electrode  kinetics 
has  been  included. 

This  model  has  been  developed  for  planar  cells  with 
cross-flow  gas  feeding  system.  It  takes  account  of  mass, 
energy  and  momentum  balances  of  the  gaseous  streams, 
together  with  the  energy  balance  of  the  solid,  considering 
the  following  chemical  and  electrochemical  reactions: 

C02  “E  y  02  T-  2e~  — >  C032-  (cathode) 
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H2  +  C032-  ->  H20  +  C02  +  2e~  (anode) 

H20  +  CO  — •>  H2  +  C02  (anode) 

The  simulation  results  are  available  in  terms  of  maps  of 
gas  and  solid  temperatures,  electrical  current,  Nernst  vol¬ 
tage,  polarisation,  internal  resistance,  pressure  drop  and 
composition  of  the  gaseous  streams  on  the  cell  plane. 

The  kinetics  is  expressed  as  a  local  electrical  resistance 
function  of  temperature,  pressure,  gas  composition  and 
phenomenological  parameters  identified  by  experimental 
data.  In  this  way,  cell  performance  is  simulated  with  a  linear 
dependence  of  cell  voltage  on  current  density,  obtaining  a 
good  agreement  between  simulation  and  experimental 
results  when  rich  reactant  feedings  are  used  [1], 

However,  experimental  data  collected  on  MCFC  stacks 
working  at  high  reactant  utilisation  factors  show  at  high 
current  densities  an  abrupt  voltage  drop  not  foreseen  by  this 
micromodel  [2].  This  trend  can  be  explained  in  terms  of 
limiting  transport  phenomena  in  the  porous  electrodes,  which 
cause  cell  voltage  decay  indicating  a  limiting  current  density. 

Aim  of  this  paper  is  to  investigate  this  limit  performance 
of  MCFC  cells  and  to  extend  to  a  wider  operating  validity 
range  the  electrode  kinetics  micromodel  accounting  also  for 
the  effects  of  the  diffusion  phenomena.  For  this  reason,  a 
number  of  specific  performance  curves  have  been  carried  out 
on  MCFC  cells  produced  by  Ansaldo  Ricerche  (ARI) 
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Nomenclature 

A,  B,  .viR,  D,  G 

parameters  in  Eq.  (2) 

c 

reactant  concentration,  mol  m-3 

E 

thermodynamic  voltage,  V 

F 

Faraday’s  constant,  C  mol-1 

J 

local  current  density,  A  m-2 

Kc 

mass  transport  coefficient,  m  s’  1 

1 

spatial  coordinate,  m 

L 

cell  length,  m 

n 

number  of  electrons 

Po2 

partial  pressure  of  oxygen,  Pa 

Q 

molar  flow-rate  for  unit  length, 
mol  s_l  m_1 

R 

electrical  resistance,  ohm 

gas  constant,  J  mol”1  K-1 

T 

temperature,  K 

u 

utilisation  factor 

V 

cell  voltage,  V 

y 

molar  fraction 

Greek  letters 

P 

parameter  in  Eq.  (2) 

YJconc 

concentration  polarisation,  V 

V 

stoichiometric  coefficient 

Superscripts 

b 

bulk 

0 

inlet 

Subscripts 

e 

electron 

lim 

limiting  value 

r 

reactant 

tot 

total 

increasing  cell  current  density  beyond  ohmic  region  (linear 
trend)  up  to  the  maximum  value  allowed  by  diffusion  limits 
without  reaching  a  zero  cell  voltage,  that  is  up  to  the  limiting 
current  density.  Tests  have  been  performed  by  separately 
varying  the  concentration  of  reactants  (C02  or  02  at  the 
cathode  and  H2  at  the  anode)  with  respect  to  a  reference  gas 
composition  to  identify  the  dependence  of  anodic  as  well  as 
cathodic  limiting  current  density  on  operating  parameters. 

The  analysis  of  these  experimental  data  requires  parti¬ 
cular  attention  because  it  is  difficult  to  isolate  the  diffusion 
effects  of  each  single  reactant  on  voltage  drop.  In  fact, 
although  the  size  of  the  cell  is  small  (about  55  cm2),  different 
parts  of  the  cell  may  work  under  diffusion  limiting  condi¬ 
tions  due  to  different  reactants.  The  merging  of  these  effects 
in  overall  voltage  and  current  measurements  complicates  the 
identifying  of  separate  anode  and  cathode  limiting  currents 
and  cases  can  be  isolated  only  when  a  sole  reactant,  under 
high  dilution  conditions,  is  limiting  on  the  entire  cell. 

In  this  paper,  we  will  present  and  discuss  our  experimental 
and  theoretical  work  and  underline  the  usefulness  of  the 


results  for  applications  to  stack  and  MCFC  plant  process 
design. 


2.  The  kinetics  model 


This  paper  proposes  an  extension  to  a  wider  operating 
condition  validity  range  of  the  electrochemical  kinetics 
relationship 

V  =  E~RJ  (1) 

discussed  in  the  previous  paper  [1], 

In  Eq.  (1)  a  linear  dependence  of  cell  voltage  V  as  a 
function  of  local  current  density  J  was  assumed  to  calculate 
the  losses  reducing  the  thermodynamic  voltage  E.  So,  the 
electrochemical  cell  performance  was  evaluated  on  the  basis 
of  a  local  electrical  resistance  R,  due  to  ohmic  and  different 
polarisation  effects.  In  particular,  R  was  written  as  the 
following  semi-empirical  relationship 

AePlT  rrr 

R  =  ^^  +  siR+DeG'T  (2) 

Po2 

where  po2  is  the  partial  pressure  of  cathodic  oxygen,  T  the 
solid  temperature,  ,viR,  /J,  B,  G,  A  and  D  are  coefficients 
evaluated  by  experimental  data  (see  also  for  numerical 
values  [1]). 

Using  these  electrochemical  kinetics  a  good  agreement 
was  obtained  for  single  cells  as  well  as  for  cells  assembled  in 
stacks  at  different  operating  conditions.  Nevertheless,  con¬ 
sidering  the  experimental  results  collected  on  cells  operating 
at  high  current  density  or  with  poor  reactant  gases,  a  good 
agreement  was  only  obtained  by  introducing  a  corrective 
coefficient  in  terms  of  a  constant  factor  penalising  perfor¬ 
mance  at  high  utilisation  factors  [2], 

The  above  consideration  has  suggested  the  need  to  inves¬ 
tigate  diffusion  phenomena  more  deeply  in  a  wider  range 
of  operating  conditions  and  on  the  basis  of  a  theoretical 
approach. 

In  particular,  the  voltage  decay  related  to  diffusion  of  each 
reactant,  well  known  as  concentration  polarisation,  can  be 
written  as  follows  [3]: 


RgT 

^1  r.conc  r,  llT 

1lrr 


(3) 


where  Jr^m  is  the  limiting  current  density,  that  is,  the 
maximum  current  reachable  when  the  limiting  reactant 
concentration  on  the  electrode  surface  approaches  zero. 

Nevertheless,  as  the  linear  contribution  of  concentration 
polarisation  has  been  already  considered  in  parameter 
validation  of  resistance  R  (Eq.  (2)),  not  to  count  this  effect 
two  times,  only  the  non-linear  part  of  Eq.  (3)  has  been  taken 
into  account: 


ir,  cone,  non- 


RJ 

nTF 


In 


(4) 
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In  this  way,  the  following  cell  voltage  expression  has  been 
used: 


V  =  E-RJ  - 


r, cone, non- 


(5) 


where  ^r,conc,non-iin  is  the  non- linear  part  of  the  voltage  decay 
related  to  diffusion  of  each  reactant. 

In  conclusion,  neglecting  O2  diffusion  effect  on  the  basis 
of  experimental  results  (see  Section  3),  the  electrochemical 
kinetics  model  proposed  is  the  following: 


f AeB/T 

V  —  E  —  |  — g - F  SiR 

P02 


DeG/TJ  J 


RJ 

nrF 


In  1  - 


J i 


H2,lim 


^H2,lim 


+  ln  1  - 


/ 


J 


Jco2  ,lim  J  Jco2  ,lim  J 


nT  =  2 


(6) 


The  limiting  current  densities  can,  on  the  basis  of 
Faraday’s  and  Fick’s  laws  and  neglecting  effects  of  counter 
diffusion  [4],  be  evaluated  as  follows: 

Jr, lim  =  nrFKacbT  (7) 


if  Cj  is  the  reactant  concentration  in  the  bulk  and  KCl  the 
mass  transport  coefficient. 

So,  in  order  to  use  the  formulation  (6),  it  is  necessary  to 
know  the  values  of  the  mass  transport  coefficients  KCr •  They 
can  be  identified  by  standard  experimental  tests  on  an  MCFC 
cell  if  correlated  to  an  entity  of  easy  evaluation  such  as  the 
utilisation  limiting  factor  Mr  iim,  which  has  been  defined  as 
the  ratio  between  the  reactant  flow  rate  consumed  when  limit 
operating  condition  extends  all  over  the  cell  and  the  reactant 
flow  rate  fed  to  the  cell. 


In  fact,  starting  from  this  definition  and  assuming  to 
operate  with  diluted  gases,  the  local  value  of  molar  fraction 
yr  can  be  written  as  a  function  of  the  local  utilisation  factor  ur 
as  follows: 


Vr=.y°(l  -  Mr) 


(8) 


Then,  the  mass  balance  of  the  limiting  reactant  flowing 
along  an  electrode  can  be  written  in  the  following  form 
(where  /  is  the  spatial  coordinate  with  origin  at  gas  flow  inlet, 
see  Fig.  1) 


(*£  =  **>  =  ***(!-*)  (9) 

with  wr  =  0  for  /  =  0;  uT  =  Mr. um  for  l  =  L. 

After  integrating  Eq.  (9),  Kcr  can  be  easily  evaluated  by 
Mr.iim  data  as  follows: 


Kcr 


g_ 

Lcr° 


ln(l  Mrlim ) 


(10) 


In  this  way,  thanks  to  a  series  of  tests  conducted  until  limit 
operation  conditions,  MrUm  and  so  KCr  have  been  experi¬ 
mentally  identified  as  discussed  in  the  following  section. 

So,  in  comparison  with  previous  approaches  [5-9],  more 
detailed  but  difficult  to  be  experimentally  validated,  the 
above  semi-empirical  voltage  decay  expression  (Eq.  (6)) 
shows  a  complete  phenomenological  description  with  a 
minimum  number  of  measurable  parameters. 

Integrating  Eq.  (6)  into  the  macroscopic  MCFC  cell 
model,  cell  behaviour  up  to  limiting  diffusion  operating 
conditions  has  been  analysed  and  a  good  agreement  bet¬ 
ween  simulated  and  measured  results  has  been  obtained 
(see  Section  4). 


Ansaldo  Ricerche  -  Lab. 
Monocella  MCFC  assemblata 


Fuel 


Fig.  1 .  Reactant  feeding  directions  in  the  MCFC  single  cell  used  for  tests. 
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3.  Experimental 

Experimental  activity  was  carried  out  on  an  MCFC  square 
planar  single  cell  of  about  55  cm2  area,  where  the  anode 
was  a  Ni/Cr  alloyed  powder,  the  cathode  a  lithiated  nickel 
oxide  structure  and  the  electrolyte  a  lithium  carbonate  and 
potassium  carbonate  mixture  [10].  The  cell  was  tested  in  the 
experimental  facility  shown  in  Fig.  2,  where  the  cell  was 
held  in  a  vessel  and  it  was  possible  to  set-up  and  control  each 
operating  parameter  [11], 

As  above  mentioned,  a  number  of  performance  curves 
were  carried  out,  increasing  cell  current  density  beyond 
ohmic  region  (linear  trend)  up  to  the  maximum  value  allowed 
by  diffusion  limits  without  reaching  a  null  cell  voltage. 

In  particular,  two  series  of  tests  were  performed  at  923  K, 
temperature  kept  constant  on  the  cell  plane  using  heating 
plates.  A  first  one  aimed  at  analysing  effect  of  reactant 
diffusion  by  varying  the  concentration  of  H2  with  respect  to 
the  reference  gas  composition  reported  in  Table  1  as  “Anode 
I”  coupled  with  a  poor  cathodic  feeding  (“Cathode  IA”)  as 
well  as  a  rich  cathodic  feeding  (“Cathode  IB”). 

In  these  tests  H2  molar  fraction  was  changed  from  0.06  up 
to  0.7  balancing  this  variation  by  reducing  N2,  as  it  is  an  inert 
component,  from  0.8  down  to  0.16  in  order  to  isolate  the  H2 
effect.  Fig.  3  shows  the  related  experimental  characteristic 
curves  obtained  at  atmospheric  pressure  using  the  “IA” 


Fig.  2.  MCFC  single  cell  test  facility. 


Table  1 

Reference  gas  compositions  and  flow  rates  (mol  s  1  X  105) 


Operating 

conditions 

co2 

h2 

h2o 

n2 

02 

Anode  I 

1.4 

6.5 

2.3 

15.8 

_ 

Cathode  IA 

14.3 

- 

- 

125.0 

15.1 

Cathode  IB 

48.0 

- 

- 

93.3 

24.8 

Anode  II 

2.4 

36.6 

3.8 

- 

- 

Cathode  II 

23.3 

- 

- 

205.4 

24.8 

composition  at  the  cathode:  the  linear  part  of  each  curve 
is  parallel  to  the  other  ones,  according  to  the  assumed 
independence  of  electrical  resistance  from  H2  concentration 
(see  Eq.  (2)).  On  the  contrary,  at  high  utilisation  factor, 
performance  decreases  disclosing  different  limiting  current 
densities  for  different  H2  concentrations. 

With  reference  to  Fig.  3,  the  current  density  values 
corresponding  to  the  minimum  voltage  achieved  along 
characteristic  curves  identify  the  best  experimental  evalua¬ 
tion  of  the  limiting  current  densities.  The  latter  are  reported 
in  Fig.  4  versus  H2  feeding  molar  fraction  (y^2).  The  figure 
shows  that  for  y^  less  than  0.15  the  dependence  of  current 
density  versus  y^  is  linear  as  expected  from  the  following 
overall  cell  balance 


7 r.lim 


Ur.lim) 


where 

Mr.lim  =  1  -  exp 


Lc°KCr\ 

Qr  J 


=  constant 


(11) 


(12) 


as  obtained  from  Eq.  (10). 

In  this  region,  then,  the  cell  seems  working  in  only  anodic 
limiting  diffusion  conditions:  the  H2  limiting  utilisation 
factor  can  so  be  experimentally  identified  as  «UKiim  =0.75 


Fig.  3.  Experimental  characteristic  curves  at  different  inlet  anodic  H2 
molar  fractions  (at  650°C  and  1  atm). 
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Fig.  4.  Measured  limit  current  density  at  different  inlet  anodic  H2  molar 
fractions  (at  650°C  and  1  atm). 


Fig.  5.  Measured  limit  current  density  at  different  inlet  cathodic  C02 
molar  fractions  (at  650°C  and  1.5  atm). 


and  the  related  average  mass  transport  coefficient  is  eval¬ 
uated  from  Eq.  (10)  equal  to  =  0.011  ms-1.  This 
value  is  very  low  with  respect  to  typical  transport  coeffi¬ 
cients  in  the  homogeneous  gas  phase  on  the  electrode 
because  it  additionally  contains  the  contribution  of  all 
transport  and  chemico-physical  phenomena  of  the  reactants 
(e.g.  diffusion,  solution,  adsorption, . . .)  in  the  complex 
structure  of  the  electrode. 

For  H2  molar  fraction  higher  than  0.15  the  curve  of  Fig.  4 
bends,  because  cathodic  diffusion  also  becomes  significant 
in  this  range.  When  C02  and  02  utilisation  factors  achieve 
respectively  0.53  and  0.24,  the  limiting  current  density 
seems  to  become  independent  of  FI2  concentration  and  to 
be  bound  only  by  cathodic  gas  diffusion.  In  order  to  verify 
this  interpretation  new  tests  have  been  performed  with  a 
richer  cathodic  composition  (“Cathode  IB”).  The  experi¬ 
mental  results  show  that  limiting  current  densities  increase 
as  expected  (Fig.  4)  and  confirm  that  the  previously  observed 
flat  behaviour  is  caused  by  cathodic  limit  conditions. 

So,  preliminary  tests  have  been  carried  out  to  allow  a 
qualitative  analysis  of  the  cathodic  diffusion  limits  follow¬ 
ing  the  same  procedure  previously  used  for  anodic  tests. 

In  this  case  the  two  reactants  (C02  and  02)  involved  in  the 
cathode  reaction  have  to  be  considered.  Operating  condi¬ 
tions  “Anode  II”  and  “Cathode  II”  (Table  1)  have  been  used 
as  reference  at  anode  and  cathode  side  respectively,  C02  and 
02  have  been  varied  parametrically  one  by  one. 

During  the  first  test,  C02  molar  fraction  has  been  varied 
from  0.02  up  to  0.15.  After  replacing  reference  conditions, 
the  second  test  has  been  carried  out  changing  02  molar 
fraction  from  0.04  up  to  0.16. 

The  results  are  reported  in  Figs.  5  and  6. 

The  curve  versus  y®02  shows  an  experimental  trend 
similar  to  the  H2  one:  a  first  linear  part  and  then  a  part 
independent  of  reactant  concentration. 


With  regard  to  02,  the  characteristic  curves  have  not  been 
carried  out  up  to  limit  operating  conditions  as  advanced  as 
for  the  previous  tests,  because  low  02  concentration  sensibly 
penalises  cell  resistance  (see  Eq.  (2))  and  it  would  be 
necessary  to  reach  a  too  low  cell  voltage. 

Fig.  6  shows  the  maximum  current  densities  achieved  in 
these  tests  plot  versus  02  molar  fraction.  The  curve  does  not 
show  any  dependence  from  02  molar  fraction  except  at  very 
low  concentrations.  Moreover,  these  concentrations  are 
never  used  during  cell  standard  operation,  also  because 
cell  feeding  is  usually  rich  in  air  in  order  to  optimise 
thermal  management. 

On  the  basis  of  this  result,  02  diffusion  effects  have  been 
considered  negligible  and  the  decay  of  curve  in  Fig.  4  has 
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0.00  0.04  0.08  0.12  0.16 

02  Inlet  Molar  Fraction 

Fig.  6.  Maximum  measured  current  density  at  different  inlet  cathodic  02 
molar  fractions  (at  650°C  and  1.5  atm). 
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been  imputed  to  C02  diffusion.  In  this  way  the  asymptotic 
trend  of  Fig.  4  has  been  used  to  identify  the  C02  utilisation 
factor  about  equal  to  0.53  and  so  to  calculate  the  average 
mass  transport  coefficient  Kc, co2  =  0.035  ms  '. 

It  is  to  be  noted  that  our  theoretical  approach  (Eq.  (11)) 
refers  to  conditions  where  the  diffusion  of  only  one  com¬ 
ponent  (H2  or  C02)  is  limiting  on  all  over  the  cell.  On  the 
contrary,  our  experimental  values  are  believed  to  refer  to 
limiting  current  conditions  where  one  or  more  components 
are  controlling  (H2  and/or  C02).  So,  among  these  experi¬ 
mental  points  the  ones  belonging  to  proportional  trends  of 
current  density  versus  reactant  concentration  are  selected 
and  compared  to  Eq.  (10)  in  order  to  individuate  the  mass 
transport  coefficients. 

4.  Model  validation 

First  of  all,  we  observed  a  systematic  discrepancy  in 
single-cell  experimentation  between  experimental  and  the¬ 
oretical  open  circuit  voltages.  In  this  case  it  has  been  taken 
into  account  within  the  model,  lowering  the  Nernst  value  by 
0.055  V  in  the  local  voltage  calculation.  This  phenomenon  is 
under  study  and  is  likely  to  be  due  to  parasitic  currents 
facilitated  by  single-cell  assembling,  in  fact  it  is  negligible  in 
stacked  cells.  On  the  other  hand,  corrections  up  to  0.2  V  are 
usual  for  other  fuel  cell  types  [12]. 

Then,  thanks  to  the  kinetics  relationship  previously  dis¬ 
cussed  and  integrated  into  the  macroscopic  MCFC  cell 
model,  a  good  agreement  has  been  obtained  between  simu¬ 
lated  and  measured  points.  Fig.  7  shows,  as  an  example,  the 
results  related  to  three  performance  curves:  at  high  >,llli 
(0.58),  when  cathodic  diffusion  effects  are  prevailing,  at 
average  v^2  (0.22),  when  both  anodic  and  cathodic  diffusion 


Fig.  7.  Measured  and  calculated  characteristic  curves  at  different  H2  molar 
fractions  (“Anode  I”  and  “Cathode  I”). 


Fig.  8.  Comparison  between  simulated  and  measured  current  densities  at 
different  operating  conditions  (indicated  by  different  symbols). 

phenomena  are  significant,  and  at  low  y^  (0.08),  when  H2  is 
the  limiting  reactant. 

Then,  the  reliability  of  the  model  has  been  verified  in  the 
overall  experimental  range  of  H2  concentration.  Fig.  8 
compares  simulated  and  experimental  current  densities  of 
the  first  test  series  in  a  parity  plot  always  showing  a 
satisfactory  agreement.  In  particular,  an  average  error  less 
than  7%  has  been  calculated,  considering  also  data  obtained 
in  diffusion  limit  conditions.  These  results  confirm  the 
reliability  of  the  theoretical  approach  of  the  model,  while 
additional  experimental  tests  are  scheduled  to  estimate  the 
mass  transport  coefficients  more  accurately  in  a  wider  range 
of  operating  conditions  (compositions,  temperature,  pres¬ 
sure,  specific  flow  rates,  . . .). 

Figs.  9  and  10  show  the  calculated  maps  on  the  cell  plane 
of  the  ratio  JIJr,iim  between  local  current  density  and  anodic 
and  cathodic  limiting  current  density,  respectively  at  the 
operating  conditions  reported  below  the  figures. 

These  operating  conditions  correspond  to  a  point  belong¬ 
ing  to  the  final  linear  trend  of  the  second  characteristic  curve 
reported  in  Fig.  6.  From  the  map  analysis,  it  can  be  observed 
that,  even  if  the  cell  is  running  under  apparently  safe 
conditions,  parts  of  the  cell  work  in  diffusion  limit  operating 
conditions  (J / —■ ►  1)-  In  addition,  on  the  cell  plane  it  is 
possible  to  detect  contemporary  an  anodic  diffusion  control 
at  the  fuel  outlet  and  a  cathodic  one  at  the  cell  corner  where 
fuel  is  fresh  and  oxidant  exhaust.  The  location  of  these 
regions  depends  on  reactant  concentration  and  local  current 
density  maps. 

Moreover,  it  is  interesting  to  observe  that  local  diffusion 
limit  conditions  can  be  reached  also  when  concentration 
polarization  values  are  significantly  lower  than  cell  voltage. 
For  example,  in  the  case  discussed  above  the  local  maximum 
r] cone on  the  cell  plane  is  only  1/5  of  cell  voltage,  but  involves 
current  density  close  to  the  limit  value. 
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Fuel  Flow  Direction  [m] 


Fig.  9.  Calculated  map  on  the  cell  plane  of  the  ratio  between  current 
density  and  anodic  limit  current  density  (Anode:  C02  1.4  x  10  5  mol  s  1, 
H2  9.1  x  10~5  mol  s_1,  H20  2.3  x  10-5  mol  s_1,  N2  13.2  x  10~5  mol  s_1; 
Cathode:  C02  14.3  X  10~5  mol  s_1,  N2  125  X  10~5  mol  s_1,  02 

15.1  X  10~5  mol  s-1;  Cell  voltage:  0.523  V). 


The  knowledge  of  J/Jr  iim  maps  is  important  in  choosing 
safe  conditions  available  for  the  entire  cell,  in  fact,  standard 
working  points  showing  local  limiting  diffusion  conditions 
have  to  be  avoided  in  order  to  preserve  cell  performance. 
This  problem  could  be  very  significant  when  stacked  cells 
are  integrated  in  power  plants  with  strong  recycles  that 
involve  the  feeding  of  diluted  gases  to  the  cells. 


Fuel  Flow  Direction  [m] 

Fig.  10.  Calculated  map  on  the  cell  plane  of  the  ratio  between  current 
density  and  cathodic  limit  current  density  (Anode:  C02  1.4  x  1 0  5 
mol  s-1,  H2  9.1  x  10~5  mol  s_1,  H20  2.3  x  10~5  mol  s_1,  N2  13.2x 
10_5mols_1;  Cathode:  C02  14.3  x  10~5  mol  s_1,  N2  125  x  10~5 
mol  s-1,  02  15.1  X  10~5  mol  s~!;  Cell  voltage:  0.523  V). 


5.  Conclusion 

In  this  paper  the  limit  performance  of  MCFC  cells  at  high 
reactant  utilisation  factors  have  been  studied.  Characteristic 
curves  beyond  the  linear  ohmic  region  have  been  obtained 
by  specific  tests;  experimental  results  have  been  analysed, 
detecting  for  reactants  the  maximum  admissible  utilisation 
factor  and  the  average  mass  transport  coefficients. 

An  original  expression  of  voltage  decay  which  considers 
the  effects  of  the  diffusion  limit  region  has  been  set-up  and 
then  integrated  into  the  MCFC  cell  model.  This  new  kinetics 
relationship  depends  on  concentration  polarisation  due  to  H2 
and  C02,  while  02  diffusion  effects  have  been  considered 
negligible. 

The  model  has  been  validated  in  a  wide  range  of  operating 
conditions,  obtaining  a  good  agreement  with  experimental 
data. 

In  particular,  when  cells  work  under  apparently  safe 
conditions,  our  model  allows  to  individuate  possible  local 
running  in  diffusion  limit  conditions  characterised  by  local 
electrical  current  close  to  the  anodic  or  cathodic  limit  value, 
even  if  concentration  polarisation  is  not  necessarily  a  large 
fraction  of  the  cell  voltage. 

So  we  underline  as  in  defining  the  global  strategies  of  cell 
control,  where  several  constraints  have  to  be  observed  at  the 
same  time  to  guarantee  safe  cell  running,  the  correct  evalua¬ 
tion  of  limiting  diffusion  conditions  also  has  to  be  taken  into 
account  together  with,  for  example,  the  estimate  of  local 
temperature  peaks  and  high  pressure  drops  [2]. 

Finally,  additional  experimental  tests  are  scheduled  to 
more  accurately  estimate  the  mass  transport  coefficients  of 
the  reactants,  which  takes  into  account  the  contribution  of  all 
transport  and  chemico-physical  phenomena  in  the  complex 
structure  of  electrodes. 
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